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Ways of improving the representation of the velocity distribution function in the solution of the 
non-linear Vlasov equation are discussed. Hermite polynomial coefficients for temperature renormal- 
ization are derived. Transformation properties of Hermite polynomials necessary to do this are dis- 
cussed. Analytic expressions for the truncation error of the temperature renormalized functions are 
compared with the computer results. 
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Introduction 

This paper deals with some mathematical aspects of the numerical solution of the non-linear 
Vlasov equation by means of eigenfunction expansions [1, 2, 3]. 1 The model studied is a one- 
dimensional, electrically neutral plasma undergoing longitudinal oscillations. Periodic boundary 
conditions are imposed on the electron density. The initial distribution function is given in dimen- 
sionless form [8] by the expression, 



/(*, v, * = 0) = (l + a cos kx)e- r2 ' 2 /V27T 

where a is the non-linearity parameter and k = 27r/L, L being the periodic length. The electric 
field is given by the Poisson equation and has the form 



E(x, t) = J dx \ f{x' , v, t)dv-x. 



The initial value of the electric field is thus E(x, t = 0) = (a/k) sin kx. The ratio (a/k) determines 
the initial field amplitude and the initial potential energy content of the system. 

In applying the eigenfunction expansion method to the numerical solution of nonlinear partial 
differential equations, it is very important to obtain as much analytic information as possible con- 
cerning expansion properties of the solution, the truncation error incurred by neglecting all but 
a finite number of terms in the expansion and the effect of truncation errors on the validity of the 
solution. Procedures of questionable validity are those in which the solution is expanded into a 
set of orthogonal polynomials without some knowledge of convergence properties of the expansion. 
When there is no proof that the solution exists on the interval over which the equation is being 
integrated, extreme caution should be exercised in evaluating numerical results. The lack of general 
methods of dealing with nonlinear partial differential equations makes one's intuition an unreliable 
guide in judging numerical results. 



1 Figures in brackets indicate literature references at the end of this paper. 
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Truncation of an expansion with time-dependent coefficients is a particularly difficult problem 
since the function that was well represented initially by a certain expansion may no longer be 
adequately represented by the same number of terms at a later time. The usual numerical tests 
employed in finite difference methods, such as conservation of certain time-independent quantities, 
may be inapplicable or not very significant. The physically observable quantities like the number 
density and macroscopic currents are determined by the first few moments of the distribution 
function of the phase space coordinates, f(x, v, t). In this sense, the first few moments of an 
eigenfunction expansion are fixed [4] since they must reproduce the observed quantities. For 
example, conservation of density in phase space and of total energy in the nonlinear Vlasov 
equation puts conditions on the zeroth, first, and second moments of the Hermite polynomial 
expansion. Of these, the zeroth is automatically conserved, while the other two do not indicate 
any trouble until long after the numerical solution has become invalid. 

Associated with the question of truncation is the problem of finding a transformation to a 
different set of eigenfunctions in the course of the calculation which would give a better representa- 
tion of the function as it evolves with time than the initial set of eigenfunctions. If such a trans- 
formation could be automatically found the truncation error would be reduced. This is the problem 
discussed below. 

In section 1 is derived a transformation of the arguments of Hermite polynomials with the aim 

of reducing the truncation error with a fixed number of terms in the expansion. The transformation 

is applied to the problem of temperature renormalization in section 2 and the expansion coefficients 

are explicitly given. 

ex. 
In section 3 error bounds are derived for the truncation and it is shown that for— ^ 1 tempera- 

k 

ture renormalization does not result in a significant reduction of truncation error. 

Comparison of analytic results with those obtained on the computer is given in appendix II. 



1. Transformation Properties of Hermite Polynomials 

This section deals with certain transformation properties of Hermite polynomials. These prop- 
erties will be used in evaluating expansion coefficients of the two-stream distribution function and 
in temperature normalization. A few necessary formulas pertaining to Hermite polynomials [5] 
will be given in eqs (1.1), (1.2), and (1.3). Hermite polynomials 2 are defined by means of the relation 



From eq (1.1) it follows that 



. rs/2= jffe^)£l. (11) 



[n/2] (—l)k n l 



where [n/2] is the greatest integer. 

From the same equation the following relation can be derived 



[n/2] n \ 

x "= 2 ww* Ion !#*-»(*)• (1-3) 



- 2 fr A!(n-2A;)! 



2 For discussion of differences in notation and generating functions see appendix I. 
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1.1. Scale Transformation 

Given a Hermite polynomial, He,\ —j= J, we would like to express it in terms of V He,„(v). 
Using eq (1.2) we obtain 

/ v \\m (-l)*n! v »-** 

\VJ £ 2 k kl(n-2k)\ '«i-«»/» (L4) 

_ 1 W (-l)^!e* n _ 2k 
e"l\2? o 2W(n-2k)\ V 

By eq (1.3) we have: 

„ / v \ 1 tefiJ (-l) fc n!e* K»-2fe)/2] (n-2A)! „ , , 

M^r^g. 2*W-2*)! g 2<l\(n-*k-2l)\ He *^ M ' (L5) 

Let k-\-l = m. We obtain 

/ v \_ 1 l^l -j (-l)Vn! 



But the summation on A: is a binomial expansion of — : (1 — eV 
This leads to the final expression 



/ V X 1 («/2] ^f(l_ e )/« 

^HV^7 = ^ S, 2»m!(n-2m)! ^—-W- d-7) 

which for e= 1 yields the required identity 

He„(^=) = He n (v), 6=1. 

1.2. Change of Origin 

Given a Hermite polynomial, He n I — -=- J, we shall express it in terms of ^ He m (v). 

Let 

f a 

v ' = — p and a = — 7= (1.8) 

Ve \Te 

From eq (1.2) we have 

[n/2] /_l\fc w ! 



[n/2] n-2k ( — l) /v /T; T 

V V tUll v 'n-2k-mi— a >\n 

A £2 k k\m\(n-2k-m)\ V 



A=0 //e = 0' 
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Summation indices in the nested double sum in eq (1.9) can be interchanged to give the following 
expression: 

u / , ,x ^ [(n ^ )/2 ^ (-I) 70 (n-m)\ , ,. rc!(-a') m 



Using eq (1.2) with ft replaced by ft — m we obtain 



,. ft!(-a') w 



He n (v'-a>) = J j He n - m (v>) m \ n _ m) y (1.10) 



///=0 



Similarly, we have 

« ft ! a' m 

He,(r'+«')-2 o »*-(«'') i^r^of (L11) 

Substituting —j= for 1/ and -7= for a' we get with the aid of eq (1.7) 

Ve Ve 

"\ VT/~ „?. »!(B-m)! 7-"=^ fc ? 2%!( n -m-2;t)! i/e "- B - 2fc(t,) - (L12) 

Let m + 2k = p 

„ (v-a\ 1 ^ 1^*1 Blog-'*(-l)g-»*(l-6)«*/ > „ , x 

«e„ — =- = — V V 2 ' i '- He»-p(v) (1.13) 

/ _ a \/>-2A- 



1 J, n!(e-l)^ tegl VVe^I, 

> /2 p ^o PHn-p)\ en ~ Av) h 2*k\{p-2k)\ 



The last sum in eq (1.13) is a Hermite polynomial of the argument — , and we obtain: 

„ (v-a\ 1 " ft!(e-l)^ /2 „ , W7 ( -a \ 

Similarly, we have 

HeJ — 1=- )= y — — He n - P (v)He p , . (1.15) 

With (1.14) and (1.15) we find, as we must, that in the limit a = they reduce to (1*7) and for e= 1 
they give eqs (1.10) and (1.11) respectively. It should be noted that imaginary values of Ve — 1 for 
e < 1 in the argument of He p are immaterial to the analysis since the product (e—l) pl2 He p 



v 



= J is always real. This can be seen by expressing He p as in eq (1.1). 

E — 1 / 
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2. Temperature Renormalization 

In the course of time the initial potential energy of a collisionless plasma is converted into 
kinetic energy by the Landau damping [6] mechanism. The velocity distribution function initially 
given by [1] 

f(v,t = 0)=^=He {) -^. (2.1) 

will change to reflect the increase in kinetic energy. The wings of the function will rise and the 
rise will be given by the set of Hermite polynomials {C n He n (v)}. It must be noted that since no 
collisions are present, the system will not be in thermal equilibrium at t ^ and consequently 
cannot be represented as it was initially, as a product of a Maxwellian and He . An exponential 
which decreases less rapidly with increasing v may however represent the velocity distribution 
function better than the set {C n He n (v)e' v212 } . We shall let 

V277V KZ " A) 

where the value of e(t) is determined by the increase in kinetic energy at the time t and C n is the 
set of coefficients associated with the new envelope function e~ r2/2€(0 . The new representation may 
be more rapidly convergent and hence require fewer terms for a given truncation error, thus reduc- 
ing the computer time, or generate a smaller truncation error for the same number of terms. 

We shall give below a transformation which will permit the use of the modified exponential 
of eq (2.2) for the envelope of the velocity distribution function. Let us assume that the velocity 
distribution function, at a given time £,.is given by the expressions 

m v j_ 

f(v, 0= ]£ £»»(*) He m (v)e *{t) 

m=0 

(2.3) 
f(v, = 2 C >MHe„{v)e mT) 



n=0 



where e(r) refers to the value of e calculated at some previous time r and e(t) reflects the current 
increase in kinetic energy. Thus we get 

\ m r oc v- v 2 v* 

C n = , ~ Y C m (t) He n (v) He m (v)e 2 ^ *(0 2 dv. (2.4) 

V27T Jlj m=o J-x 

Let e(T)e(Q+e(Q-€(T) 

let 

vvy = v'. 

Then 

C„(0 = ~^=~ 4= 2 C "< f" He " (-7=) He " (-7=) e-'-'^dv' (2.5) 

V27m! Vy & J-> Wy/ Wy/ 

and conversely we have 

C m (t) -L- .4= V C„(0 l X He n (-7=) He m (4=-) e-»"'W. (2.6) 
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In the above equation 



6(t)6(T)+€(T)-6(t) 

P e(t)«(r) 



and 

vVp=v'. 

Let us assume that the C m are known. Using eq (1.7) to express He n (^=) m term s of He n {v') we 
obtain with the use of orthogonality conditions the following expression for C n (t): 

I m [w/2] [n/2] i m!(l— y) A ' +/ 



" (t)_ V^ m ?o £ S Cm(0 r ( "' + " )/2 2^'/t!(m-2yt)!/!( ra -2/)! 8 "- 2(,ffl - 2fc - (2 " 7) 



The orthogonality condition requires that 2k = m — n + 2l and eq (2.7) reduces to 

1 » [n/2] m l(l _y)[(m-»)/2] + 2/ 

C " (t)= V^,„? Cm(f) ,? r ( " , + " )/2 2 < »'-" +4 '» 2 ((m-» + 2/)/2)!((n-2/)!) 2 /! 



(2.8) 



with the negative factorial being denned as zero. The double sum in eq (2.8) can be written in a 
different form in which no negative factorials occur: 

1 [»/2] M m \ (1 _y)[(m-w)/2] + 2Z 

Cn{t) = V^ ,? «,?.„ Cm(0 y (,M + ll)/2 2<' M - w + 4 ')/ 2 ((/ii-n-h2/)/2)!((n-2/)!) 2 /! (2 ' 9) 



It is very difficult to make any comparison between C n and C m without the knowledge of the actual 
values of C m . Although a bound can be given for C n in terms of C under the assumption that the 
C m are monotonically decreasing with increasing m, this bound is of no practical value because 
it is too generous. 

We shall now turn to the simpler problem of calculating the coefficients in the expansion 



In this case we have 



e-^=^CAW,e-^ 2 . (2.10) 

n = 

c(0 = l; €(t) = c; y = - 
C =l; C m = 0, m#0. (2.11) 



Thus the result in eq (2.9) applies to this case for the values of y and C given above. Because the 
functions in eq (2.10) are symmetric in v, only Hermite polynomials with even indices will be 
present in the expansion. Substituting the values of y and C m into eq (2.9) we obtain: 
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3. Error Bounds for the Expansion 

We can obtain an error bound on the truncation of the expansion after the first TV terms. Using 
the asymptotic expression for Hermite polynomials [7] (see appendix I, eq (1.6)) we obtain: 

C«He tq (v)~ 2 , r( ;V ^T (— J • (3-D 

Using Stirling's formula eq (3.1) reduces to 

C 2Q lHe- iq (v) ~ Jl e^ (e ~* )g . (3.2) 

The truncation error is given by the following expression: 



^'AvrV+i+j)" 1 



Z, o Vt7^ (AT+l+p)^- 
We know that 



(3.3) 



(/V+l+p) 1 ' 2 (/V-hl) 1 / 2 
for all values of p in eq (3.3) because e > 1. 



(3.4) 



R™ < Vf ^ 7^i n,l_, (3.5) 



(6-l)- v+1 

f 7r~ (7V-hl) 1 / 2 (2-e) 



where the summation over p in eq (3.3) has been performed using (3.4). 

Let us now estimate a representative value of e that may occur in the solution of non-linear 
Ylasov equation. From the introduction we recall that the initial distribution function is given by 
the expression 

/(*, v, t = 0) = (1-ha cos kx)e'- v ' 2l2 lV27T (3.6) 

where k = 27r/L and L is the length of the system. The initial kinetic energy content of the system 
is given by the integral 

K= \ dx v 2 f(x, v, t = 0)dv = L. (3.7) 

The initial potential energy content is given by 

fL CL a 2 

P= E 2 (x)dx = — sin 2 kx dx (3.8) 

In eq (3.8) the electric field was obtained from eq (3.7) with the aid of Poisson's equation in dimen- 
sionless form. Integrating eq (3.8) we obtain the initial potential energy content 

P=$L (3.9) 
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20 40 60 80 100 120 140 

(N+-I) 

N 

Figure 1. Truncation error bound B N+/ versus (N + l) /or the function e~ v2/2€ = V C 2 n|He 2 n(v)/br v=2.0, 

6=1.31 and e=1.80. 
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Choosing a = 0.2 and k = 0.2 we obtain 

P = 0.5L (3.10) 

After the electric field has been completely damped out and the initial potential energy of the field 
converted to kinetic energy, e can be calculated from the equation 



JLf 



1.5 = — t= t?e- v2l2t dv (3.11) 

For the values of a and k given above we obtain the value of e 

e=1.31037. (3.12) 

Figure 1 is a plot of the bound Bn+i on the truncation error, Rn+i, as a function of the number 
of terms, N, for v = 2.0, e = 1.31. Figure 2 is a comparison of the number of terms necessary to repre- 
sent the function given in eq (2.10) with a truncation error Rs+\ < 10~ 5 T obtained analytically and 
numerically (see appendix II). The value of velocity chosen for the graph is v = 2. As can be seen 
by comparing the analytic and the computer curves, the values of e close to 1 are in error because 
the asymptotic form of Hermite polynomials does not apply. One thing that can be concluded from 
the graphs, however, is that the truncation error cannot be reduced appreciably by reexpressing 
the solution of the Vlasov equation by a function of the type given in eq (2.10) for the ratio a/k ^ 1. 
For ratios a/k > 1, the potential energy content of the system will lead to a larger increase in kinetic 
energy as the electric field is damped and a transformation of the type discussed would result in 
a better representation of the solution. 



4. Appendix I 

There are two different kinds of Hermite polynomials in the literature. One of these is based 
on the generating function 

e ^-, = ^tLMH (L1) 

n=0 "" 

and the other, used in this paper, is usually designated in the literature by He [6]. It is based 
on the generating function 

^ He„(x)t" 

n=0 n ' 

The asymptotic expression in the literature [7] is given in the form 

V(in+l) exp(-i*; 2 ) //„(*)= T(n+1) [cos(N l ' 2 x-inir) +0(rc- 1/2 ) ] (1.3) 

where N=2n+1. 

There is a simple transformation that connects the two polynomials. Let 

t' x 

^7T and * = vr' d.4) 
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In eq (1.1) we obtain: 



e 2, 2 ^n\ A n! {Li)} 



The last equality follows from eq. (1.2). 

Using the transformed variable of eq (1.5), the asymptotic form for He,, (a;) is 



/M*)= eX2,4r(n+1) 



2»/ 2 r(i7i+i) 

where again N = 2n + 1 . 



cos (^ 1/2 ^-f ^) + °("~ 1/2 ) 



(1.6) 



5. Appendix II 

A computer program was written in Fortran to calculate the terms of the expansion of/(v, t) 
of eq (2.10). The He 2n (v) and C 2n were calculated recursively. Scaled polynomials, He 2n , were used 
to stay within the range of the computer in single precision. 

He 2n (v) =—==He 2 n(v) (II.l) 

V(2n)! 



C 2 „ = V(2/i)!C 2 «. 

A partial sum, P M , of the first n terms was formed and more terms were added until the relative 
error ratio test of eq (II. 2) was satisfied for one hundred successive terms. 

| R.E. (P. ) | = |P ' , :i"T ,P " 1 < 10- 5 . (II.2) 

Using eqs (3.2), (3.5), and (II. 2) the three term inequality relating Rs+i to the truncation error bound 
107 5 is: 

Rn + i < ; s r" < l0 ~ 5T (113) 

V / ]V+T(2-€) 

VJV+1 (2-6) 
where /is the function of eq (2.10). The parameter values investigated were 

e=1.10, 1.20, 1.40, 1.60, 1.80 

z; = 0.00 to 6.00 in intervals of 0.25. 

Figure 2 contains a computer curve for i? = 2.0 and a comparison curve based on eq (II. 3). 

It should be noted that a uniform bound on the truncation error can be obtained by allowing 
v to assume its largest value. 
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FIGURE 2. N versus e for v = 2.0 for temperature renormalization curve. 
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